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1.  Introduction 


The  extreme  values  of  a  stochastic  process  rather  than  its  usual  values  are 
often  of  paramount  interest.  In  a  manufacturing  plant.  Lor  example,  a  typical 
concern  is  the  queue  length  of  parts  waiting  to  be  processed  at  a  work  station. 
Small  to  moderate  values  of  the  queue  may  indicate  that  the  system  is  operating 
successfully  and  fluctuations  in  the  queue  are  unimportant.  On  the  other  hand, 
large  queues  may  call  for  extraordinary  measures  such  as  the  allocation  of 
auxiliary  storage  space,  employee  overtime,  or  rescheduling  of  production.  A 
natural  question  is:  What  is  the  probability  that  the  queue  will  exceed  a 
specific  critical  value  in  a  certain  time  period?  Such  questions  are  the  focus 
of  this  study.  The  following  paragraphs  describe  the  gist  of  our  results. 

Consider  an  M/ 1/1  queueing  system  in  which  customers  arrive  at  a  single 
server  according  to  a  Poisson  process  with  rate  \  and  the  independent  identi¬ 
cally  distributed  service  times  have  a  mean  y  \  The  queue-length 

process  is  positive  recurrent,  null  recurrent,  or  transient  according  to 

whether  the  traffic  intensity  p  =  A/y  is  below,  equator  greater  than  1, 

respectively.  Let  denote  the  maximum  queue  length  up  to  the  n^  time  the 

system  becomes  empty.  We  say  that  M  has  the  limit  distribution  F  if  there 

are  norming  constants  a  ,b  >0  such  that 

n  n 

lim  P ( (M  -a  )/b  <  x)  =  F(x), 

n  n  n  — 

n-K» 

for  each  continuity  point  of  F.  If  there  are  no  such  F,  a  ,b  ,  we  sav  that  M 
y  n  n  n 

does  not  have  a  limit  distribution. 

From  the  classical  extreme  value  theory,  as  in  Galambos  (1978)  or  Leadbetter 
et  al.  (1983),  we  have  the  following  results: 


(1.1)  If  P  =  1,  then  lim  P(M  /n  _<  x)  =  e  ,  x  >_  0. 

n-+°° 

(1.2)  If  p  <  1,  then  M  does  not  have  a  limit  distribution. 

n 

Cohen  (1969)  proved  (1.1)  and  an  analogue  for  GI/M/1  queues  as  well. 

Related  results  appear  in  Iglehart  (1972)  for  GI/G/1  queues  and  in  Serfozo 

(1986)  for  birth  and  death  processes.  Other  insights  on  extremes  of  queues 

appear  in  Heathcote  (1965)  and  Heyde  (1971).  Anderson  (1970)  observed  the 

anomaly  (1.2)  and  gave  bounds  on  the  distribution  of  M^.  The  result  (1.2)  is  due 

to  the  property  that  if  is  the  maximum  of  n  independent  integer-valued  random 

variables  with  distribution  F,  then  M  does  not  have  a  limit  distribution  when 

n 

[ F (m)  -  F (m-1) ] / [ 1  -  F(m-l)]  does  not  converge  to  1;  see  for  instance  Theorem  1.7.13 
of  Leadbetter  et  al.  (1983).  The  results  (1.1)  and  (1.2)  say  that  the  maxima 
of  null  recurrent  queues  have  a  limit  while  the  maxima  of  positive  recurrent 
queues  do  not.  This  is  the  reverse  of  what  one  would  anticipate. 

In  this  paper,  we  show  that  maxima  of  positive  recurrent  queues  are  more 
well  behaved  than  (1.2)  suggests.  The  key  idea  is  that  by  allowing  the  basic 
parameters  of  the  queueing  process  to  vary  with  n  in  a  certain  manner,  then  their 
maxima  do  indeed  have  the  three  limiting  distributions  shown  in  (2.1).  The 
first  and  third  distributions  are  the  classical  Frechet  and  Gumbel  extreme-value 
distributions,  and  the  other  distribution  is  a  new  one.  We  have  shown  in  Serfozo 
(1986),  that  the  maxima  of  positive  recurrent  birth  and  death  processes  and  M/M/s 
queues  also  have  these  limit  distributions.  The  present  analysis  differs  in  that 
it  involves  subtle  arguments  dealing  with  the  convergence  of  integrals  of  complex¬ 
valued  functions.  In  Section  2  we  discuss  M/G/l  queues  and  in  Section  3  we 
discuss  GI/M/1  queues. 
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MW 
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V.'vVv 


2.  Results  for  M/G/l  Queues 


Property  (1.1)  encouraged  us  to  investigate  the  maximum  M  of  positive 

n 


recurrent  M/G/l  queues  when  the  arrival  rate  A  and  service  distribution  G  vary 
with  n  such  that  the  traffic  intensity  is  nearly  1.  This  is  the  basis  of  the 
following  results. 


Consider  a  sequence  of  positive  recurrent  M/G/l  queues  indexed  by  n  =  1,2, 


th 


where  the  n  process  has  arrival  rate  A  and  service  distribution  G  with  mean 

n  n 


y  \  For  the  n^  process,  let  M^  denote  the  maximum  queue  length  up  to  the  n^ 


time  the  system  becomes  empty.  With  no  loss  in  generality,  assume  that  the 
system  is  empty  at  time  0.  We  shall  study  the  convergence  in  distribution  of 


(M  -a  )/b  as  the  traffic  intensity  p  =  A  /y  tends  to  1.  Although,  our 
nnn  nnn  & 


limiting  results  are  for  large  p^,  the  limits  yield  good  approximations  for  the 


distribution  of  M  for  small  as  well  as  large  p  . 

n  n 


We  assume  that  the  Laplace  transform 


Y„U>  - 


e-ztdC  (t). 


for  z  complex,  is  analytic  on  the  region  Rez  >  -x  ,  where  -°°  £  -x  <0,  and  that 


Y  (x)-  -<•  oo  for  real  x  1  -x  .  This  assumption  is  equivalent  to  y  (x)  being  finite 
n  n  n 


in  some  real  interval  |x  <  e  and 


-x  =  inf  {x  <  0:  y  (x)  <  °°}. 
n  —  n 


One  consequence  is  that  all  moments  of  G  exist.  We  define 

n 


0  (z)  =  y  (A  ( 1— z ) )  -  z,  for  Rez  <  1  +  x  /A  . 
nnn  ’  n  n 

Lemma  2.2  says  that  4>  has  exactly  one  zero  on  the  real  interval  (1,  1+x  /A  ). 

n  n  n 

We  let  r^  denote  this  zero.  Let  G^  be  another  such  distribution  with  Xq,  <£q 


and  Tq  defined  similarly. 


We  will  show,  under  some  natural  conditions,  that  if  n(r  -1)  c,  then  M 

n  n 

has  the  limit  distribution  F  ,  where 

c 


(2.1) 


Fn(x)  =  exp(-x  ) 


when  c  =  0 


F^(x)  =  exp[-c/(e  -1)]  x>0,  when  0  <  c  <  00 


Fjx)  =  exp  (-e  ) 


when  c  =  00 . 


In  particular,  we  will  show  that 


(2.2) 


where 


(2.3) 


lim  P((M  -a  )/b  <  x)  =  F  (x)  xeR, 

_ _  n  n  n  —  c 

n-«o 


b  =  {3 
n  n 


a  =  xa 
n  In 

[<*  +  6n  108  11 

'  n  n  n 


c  =  0 


0  <  c  <  00 


a  =  3  log[(l-p  )/4>*  (r  )]  3  -  1/logr  . 

nn  nnn  n  n 


THEOREM  2.1.  Suppose  that 


(2.4) 


lim  G  (x/X  )  =  Gn (x) 
n  n  u 

n-K» 


for  each  G-  continuity  point,  r_  =  1,  lim  a  ^  -00,  and  when  n(r  -  1)  -*■  0 
u  rr^30  n 


(2.5) 


lim  (1-p  )/[n(r  -l)(l+e)  ]  =  0  for  each  e  >  0. 
n  n 

n-*» 


Then  the  possible  limit  distributions  for  M  are  F  ,  0  <  c  <  °°.  The  M  has  limit 

n  c  —  —  n 

distribution  F^  if  and  only  if  ntr^-l)  -►  c.  Appropriate  norming  constants  for 
the  convergence  (2.2)  are  given  by  (2.3). 


REMARKS  2.2  (i)  The  distributions  F  ,  0  <  c  <  00  are  of  distinct  type:  F  and  F  . 

c  —  —  c  e 

are  of  the  same  type  if  and  only  if  c=c' . 

(ii)  For  any  positive  recurrent  M/G/l  queue,  without  the  artificial 

dependence  on  n,  Theorem  2.1  yields  the  approximation 


(2.6) 


P (Mnlogr  -  log[  (1-p) /<}>'  (r)  ]  <_  x)  =  exp[-n(r-l)/ (eX-l)  ] ,  x  _>0. 


There  are  analogous  approximations  of  the  distribution  of  by  F^  and  F^  when 
n(r-l)  is  very  small  or  very  large.  However  (2.6)  is  better  since  the  actual 
value  of  n(r-l)  appears  on  the  right  whereas  it  doesn't  for  the  F^  and  F^ 
approximations,  which  represent  degenerate  cases.  Furthermore,  since  n(r-l) 
appears  on  the  right,  (2.6)  is  accurate  when  r  (or  p)  is  not  near  1.  For  the 
M/M/1  queue,  we  found  that  the  difference  between  the  two  sides  of  (2.6)  is 
below  0.018  when  n  >  15  and  below  0.01  when  n  >  20  for  any  p  in  (0,1). 


(iii)  Assumption  (2.4)  says  that  G^C'/A^)  converges  weakly  to  G^ . 

This  implies  that  y  (A  z)  -*■  y_(z)  and  <p  (z)  -*■  <t>„(z)  uniformly  on  compact  sets. 

n  n  u  n  U 

We  use  the  latter  and  rn  =  1  to  ensure  that:  (a)  r  1,  and  (b)  there  is  an 

u  n 


n  >  1  and  N  such  that  4>  (z)  for  n  >  N  has  no  zeros  on  the  annulus  1  <  Izl  <  n 

n  —  11 

other  than  at  r  .  Theorem  2.1  is  true  with  (2.4)  and  r^  =  1  replaced  bv 
conditions  (a)  and  (b) . 


(iv)  We  use  assumption  lim  a  ^  - 00  along  with  r  ■+  1  to  ensure  that 

n-**  n  n 

+  b^x  00 .  Assumption  (2.5)  is  only  for  the  degenerate  case  ntr^-l)  -*■  0. 

These  two  assumptions  are  automatically  satisfied  in  natural  cases  such  as  those  in 


Corollary  2.4. 

To  prove  Theorem  2.1,  we  use  the  following  properties  of  <j>  ;  we  suppress 

n 

the  n  here  for  convenience. 


mm 


- 


PM! 


_vwv»> 


%*1 
I] 

,**! 

i*2| 


LEMMA  2.3.  (i)  The  function  <J>(x),  for  real  x,  is  strictly  convex  on  the  interval 
(0,  1+x/A)  and 

(2.7)  <t>(x)  00  as  x  f  1  +  x/A. 

Furthermore,  4>(x)  has  exactly  two  zeros  on  this  interval  at  1  and  at  some  r  >  1. 

(ii)  The  only  zeros  of  <f(z)  on  the  circles  |z|  =1  and  |z|  =  r  are  at  1  and  r, 
respectively,  and  these  zeros  are  of  order  1. 

(iii)  $(z)  has  no  zeros  on  the  annulus  1  <  lz|  <  r. 


Proof .  (i)  The  second  derivative  of  <t>(x)  exists  and  is  positive,  and  so  <f> 
is  strictly  convex.  Statement  (2.7)  is  obvious  when  x  <  °°,  and,  when  x  =  °°, 
it  follows  by  the  strict  convexity  of  <p  and 


(2.8) 


lim_  <f>'(x)  =  lim  tyet^dG(t)  -  1 
x+l+x/A  y~+°° 


Clearly  c£(l)  =  0  and  so  the  strict  convexity  and  (2.7)  ensure  that  there  is 
exactly  one  other  zero  of  (p  at  some  point  r  in  the  interval  (1,  1+x/A)  . 

(ii)  Since  y(A(l-x))  is  strictly  increasing  and  convex  with  fixed  points  at  x  = 
and  x  =  r  (the  zeros  of  <J>(x)),  then 


y(A(l-x))  <  1  when  0  <  x  <  1,  and 

y(A(l-x))  <  r  when  0  <  x  <  r. 


Then  on  the  circle  |z|  =  1,  the  only  zero  of  <)>  is  at  1,  since 


|<t(z)  |  _>  |  z|  -  |y(A(1-z))  | 

>  1  -  Y(Ml-Rez))  >  0  when  z  ^  1  (Rez  <  1) 


Similarly,  on  the  circle  |z|  =  r,  the  only  zero  of  <f>  is  at  r,  since 


| ( z )  |  >  r  -  Y(A(1-Rez))  >  0  when  z  ^  r  (Rez  <  r)  . 


The  zeros  of  $  at  1  and  r  are  of  order  one  since  $'(1)  =  p-1  <  0,  and  (r)  >  0  by 
strict  convexity. 


•**/  *VV,«»  '  J*'  ^  »'<'V4  ,""S> 


■y.' 


7 


(iii)  By  Rouche's  theorem  and 

|  4>  ( z)  -  z|  =  |y(X(1-z))  |  <_  y  ( X  ( 1-  |  z  [ ))  =  4>(|z|)  +  |  z  |  <  |  z  |  when  1  <  |z|  <  r, 

we  know  that  $(z)  and  z  have  the  same  number  of  zeros  in  1  <  |z|  <  r.  Thus  <}> 
has  no  zeros  there. 


Proof  of  Theorem  2.1.  We  shall  first  prove  that  n(r  -1)  -*  c  implies  (2.2). 
-  n 

th 

Associated  with  the  n  M/G/l  process,  let  denote  the  distribution  of  the 
maximum  queue  length  in  one  busy  period  (the  period  between  two  successive 
entrances  to  state  0).  From  Takacs  (1965)  and  Cohen  (1969),  we  know  that 

(2.9)  F  (m)  =  1  -  (2iri)  1  (l-z)z  m4>  (z)  1dz/(2’TTi)  1  z  m<{>  (z)  "^dz, 

n  J  .  n  J  .  n 

C  c1 

where  is  any  circle  in  the  z-plane  with  center  at  the  origin  and  with  radius 
less  than  1,  and  the  complex  integrals  are  over  C1  in  the  counter-clockwise 
direction. 

Because  the  queueing  process  regenerates  each  time  the  system  becomes 
empty,  then  is  the  maximum  of  n  independent  random  variables  with  distribu¬ 
tion  F^.  Consequently, 

P((M  -a  )/b  <  x)  =  F  (a  +  b  x)n. 

n  n  n  —  n  n  n 

It  is  well  known  that  (1  +  n^)0  **  e^  for  any  real  numbers  satisfying 
nHn  -*•  H  where  -°°  <  n  <  00 .  Then  to  prove  (2.2),  it  suffices  to  show  that 


when  c  =  0 


(2.10)  lira  n[l  -  F  (a  +  b  x) ]  =  c/(e  -1) 

n  n  n 

n-*°° 

-x 

e 


0  <  c  <  oo 


To  this  end,  consider  the  integrals  in  (2.9).  We  shall  express  them  in 

•.  r-orc  convenient  form.  Choose  circles  and  C  in  the  z-plane  with  centers 

n  n 


'A‘-' 


'j  *>*.«  V  v  v  .  ..  .■>  .'-V-  .V  .V  .v  .% 

1  .*  ■  f  V  *  J*  *  **  **■  •'j  J*  •  /*'  «"■***■*"*  **1  J  m**1’  „**•  j.*''  kV  i,-"  '  "*  "'*■  ^ 


at  the  origin  and  radii  less  than  1  and  greater  than  r  ,  respectively,  such 

that  d>  has  no  zeros  other  than  1  and  r  on  the  annulus  between  and  C  . 
n  n  n  n 

This  is  possible  because  of  the  nature  of  0^;  see  Lemma  2.3.  For  convenience, 
we  temporarily  suppress  the  subscript  n  on  (f)^,  C^,  C  ,  r^,  p^. 

From  Lemma  2.3,  we  know  that  the  function  z  m$(z)  ^  has  poles  of  order 
one  at  1  and  r  and  is  analytic  at  all  other  z  in  the  annulus  between  c''  and  C. 
The  residues  of  this  function  at  these  poles  are 


lim  (z-l)z  m<f>(z)_1  =  1/<J>'  (1)  =  -1/O.-0) 
z~+l 

lim  (z-r)z  m4) ( z)  ^  =  1/  (rmcf> ’  ( r ) )  . 
z-+r 

Then  by  the  residue  theorem  for  complex  integration,  it  follows  that 


(2.11) 


(2TTi) 


■  _ 

z  m<Kz)  dz 
JC 


(2-rri) 


i  z  m$(z)  1dz 
!C 


=  -1/ (1-p)  +  l/(rV(r)) 


Similarly , 

(2.12) 


(2-rri) 


-1 


(l-z)z  m<j>(z)  ''’dz  -  (2iri)  ^  ^ 


(l-z)z  mcJ> ( z )  ^dz 


(1-r)  /  (rm4> '  (r)). 


where  the  last  term  is  the  residue  of  the  integrand  of  these  integrals  at  r; 
the  residue  at  1  is  0. 

Combining  (2.9),  (2.11),  (2.12)  and  letting  m^  denote  the  integer  p'^t 

of  a  +  b  x,  we  have 
n  n 


m  -1 

n[l-F  (a  +b  x)  ]  =  [n(l-r)/(r  n$'(r))  -  n(2TTi) 
n  n  n 


-m  , 

(l-z)z  °'Kz)  dzi 


/[-!/( 1-0)  + 


^  —  1  —  1 
r  V(r)  -  (2-i) 


z  n ;  (  z )  ^  d  7.  ] . 


'  J*  f  r  1*  /  rn“  .  ,*  /  ^  _  *  m  "  ‘  . 


9 


That  is, 


(2.13)  n [ 1-F  (a  +b  x) ]  =  (1+6  )/{(?  -1 ) / [ n ( r  -1)]  +  e  } 
n  n  n  n  n  n  n 


where,  returning  the  subscript  n  to  '^n  »  rn  *pn  *  Cn » 


5  -  r  U  )/(l-Pn) 

n  n  n  n  n 


5  =  t  (1-P  )  [n(r  -1)1  1  ( 2TTi )  1 

n  n  n  n 


(l-z)z  n$  (z)  dz 
n 


1  -i  ’  “in  . 

e  =  C  (1-p  )[n(r  -1)]  (2Tii)  z  (z)  dz, 

n  n  n  n  n 

n 

We  first  note  that  n(r  -  1)  -*■  c  implies 

n 


c  =  0 


(2.14)  lim  (l,  -1)  /  [  n(r  -1)]  =  (e‘-l)/c 

n  n 

n_>00 

e 


0  <  c  <r  oo 


This  follows  since 


(2.15) 


=  exp([a  +b  x  +  0(1)]  log  r  -  ot  /B  } 
n  n  n  n  n  n 


1  +  xn(r  -1)  +  o(l) 
n 

x+o ( 1 ) 


x+o ( 1 ) 


0  '  c  <- 


n (r  - 1 ) e 
I  n 


For  case  c  -  0,  use  eU  =  l+u+o(u)  as  u  -  0  and  (v  -  1)  ^  log  v  1  as  v  4-  1 . 

Ue  now  show  that  *  and  ?;  converge  to  0.  Choose  r|  in  the  interval 
n  n 

(1,  1+Xq)  such  that  t0  has  no  zeros  on  the  annulus  1  \z\  •  n.  This  is 

possible  due  to  the  nature  of  *  Since  2  (z)  •*  yM(z)  uniformly  on  compact 
sets,  there  is  a  positive  integer  N  such  that  }  (z) ,  for  n  N,  has  no  zeros 


z  ■  r  other  than  at  r  .  Therefore,  we  can  take  each  C  ,  for  n  >  N, 

n  n  — 


to  be  the  circle  with  radius  r.  Letting  =  sup  ■  l^izl 


z !  =  n  1  and 


•  c-..  -  .  .  - 

"k  -V  ~r 


using  (2.15)  we  have 


6' 

n 


<  ?  (1-p  )[n(r  -1)]  1  (2tt)  1 
—  n  n  n 


-m 


1— z I  I  z I  n|4>  (z)  I  1dz 


I  z  |  =n 


nx+a  +0(1) 

0 (1) B  (1-p  ) / [n(r  -l)n  n  ]  c  =  0 

n  n  n 

m 

o(i)b  (i-p  )/n  n  o  <  c  <  °°. 

n  n  — 


Applying  B^  -*■  Bq  <  °°,  assumption  (2.5),  and  m^  -*■  00  to  these  expressions,  it 

follows  that  6^  ->  0.  A  similar  argument  yields  -*■  0. 

Using  these  limit  statements  for  5^,  6^  and  in  (2.13)  implies  (2.10), 

which  in  turn  yields  (2.2).  Thus,  we  have  shown  that  n(rn~l)  ■*  c  implies  (2.2) 

To  prove  the  converse,  suppose  has  the  limit  distribution  F  .  Let 

n'(rn,-l)  be  any  convergent  subsequence  of  n(rn~l)  on  the  compact  set  [0,°°] 

and  let  c'  =  lim  n' (r  ,-1).  Then,  from  what  we  iust  proved,  M  ,  has  the  limit 

k^o  n  n 

distribution  F  , .  But  M  , ,  as  a  subsequence  of  M  ,  also  has  the  limit  distri- 
c  n  n 

bution  F  .  Then  by  Khinchine's  Theorem  (see  for  instance  Theorem  1.2.3  of 

c 

Leadbetter  et  al.  (1983)),  the  distributions  F  ,  and  F  are  of  the  same  type. 

c  c 

Consequently  c'=c.  Thus,  any  convergent  subsequence  of  n(r  -1)  must  converge 

n 

to  c,  and  hence  n(rn~l)  -*•  c.  A  similar  argument  shows  that  F  ,  0  c  <  00 ,  are 
the  only  possible  limit  distributions  of  M^. 

Two  examples  of  Theorem  2.1  are  as  follows. 


COROLLARY  2.4.  The  assertions  of  Theorem  2.1  are  true  under  the  single 

hypothesis  p  -*-1  when  (i)  the  service  times  are  constant  or  (ii)  each  G 
n  n 

is  a  gamma  distribution  with  scale  parameter  n  and  order  k. 

Proof .  It  suffices  to  verify  the  hypotheses  of  Theorem  2.1  for  the  two  cases. 
First,  suppose  the  service  times  are  constant.  Then  y^(z)  =  exp(-z/u  )  and 


mm 
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4>  (z)  =  exp[-p  (l-z)]-z  for  each  complex  z  (x  =  °°) .  The  zero  r  of  on 
n  n  n  n  n 

(l,00)  is  given  by  (r  -1)  ^log  r  =  p  .  Assumption  (2.4)  is  satisfied  since 

n  n  n 

p  -*•  1  implies  that  y  (X  z)  =  exp(-p  z)  -*■  Y„(z)  =  exp(-z).  Clearly  r  -*•  r 
n  nn  n  0  n  Q 

Now,  using  4>'(z)  =  p  [ 4>  (z)  +  z]  -  1  and  $'(r  )  =  p  r  -1,  we  can  write 
n  n  n  n  n  n  n 

a  =  {log[ (r  -1)  -  log  r  ]  -  log[r  logr  -  (r  -l)]}/logr  . 
n  n  n  n  n  n  n 

Then  by  four  applications  of  L'Hospital's  rule  and  r^  1,  it  follows  that 
-*■  -1.  Also,  assumption  (2.5)  is  satisfied,  since  an  application  of 
L'Hospital's  rule  gives 


(1-p  ) / (r  -1)  =  [(r  -1)  -  log  r  ] / (r  -1)  ■+  1/2. 

n  n  n  n  n 


Next,  assume  that  is  a  gamma  distribution  with  Laplace  transform 


y  (z)  =  (1  +  n ^z)  k, 
n  n 


Rez  >  -x  =  -r\  . 
—  n  n 


Then  0  (z)  =  [1+(P  /k)(l-z)]“k  -  z,  Rez  <  1  +  n  /X  and  its  zero  r  is  one  of 
n  n  n  n  n 

the  k+1  solutions  of  1-r  [ 1— ( X  /r\  ) (r  -1)]  =0.  Assumption  (2.4)  is  satisfied 

n  n  n  n 

-1 

since  p  =kX  /n  ->1  implies  that  y  (X  z)  ■*  Yn(z)  =  (1+k  z)  .  Clearly 
n  n  n  n  n  v 

r  -*■  rQ  and  rQ-l  since  <J>g(l)  =  Pg-1  =  0.  Using 

1+k'1 

(J)'(z)  =  p  ($  (z)+z)/[l  +  p  /k)(l  -  z)  ]  -  1  and  <j>'(r  )  =  p  r1  -  ] 
nnn  n  nnnn 

we  can  write 

a  =  { log[  (r  -1)  -  k(l-r_1//k)  ]  -  log[k(r1+1/,k-r)  -  (r-l)]}/log  r  . 
n  n  n 

From  five  applications  of  L'Hopital's  rule,  it  follows  that  ot^  +  -1.  Finally, 
assumption  (2.5)  is  satisfied  since  an  application  of  L'Hospital's  rule  gives 


(1-p  ) / (r  -1)  =  l(r  -1)  -  k(l-r  _1/k)]  /  (r  -l)2 
n  n  n  n  n 


(1+1/k) /2 . 
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1 

as 


3.  Results  for  GI/M/1  Queues 

In  this  section,  we  show  that  the  limiting  behavior  of  extreme  values 
of  GI/M/1  queues  is  analagous  to  that  for  M/G/l  queues.  Consider  a  sequence 
of  positive  recurrent  GI/M/1  queueing  systems  indexed  by  n  r  1,2,...  For  the 
n*1^  system,  let  G^  denote  the  distribution  of  the  times  between  arrivals  and 
let  u  1  denote  the  mean  of  the  exponential  service  times.  Assuming  the 
system  is  empty  at  time  0,  let  M^  denote  the  maximum  queue  length  up  to  the 
nC^  time  the  system  becomes  empty. 

Suppose  that  the  Laplace  transform  y  (z)  of  G  has  the  same  form  as  in 

n  n 


Section  2,  and  let 


ip  (z)  =  y  (y  (1  -  z)>  -  z  for  Rez  <  1  +  x  /y  . 
n  n  n  n  n 


An  easy  check  shows  that  ^  has  the  same  properties  as  (J)^  in  Lemma  2.2  with 

in  place  of  and  rn  and  1  reversed.  In  particular,  ^(x)  has  exactly  two 

zeros  on  (0,  1+x  /u  )  at  1  and  at  some  r  <  1,  and  W  (z)  has  no  zeros  in  the 
’  n  n  n  ’  rn 

annulus  r  <  Izl  <  1. 
n  '  ’ 

For  the  following  result,  we  use  the  notation  in  (2.2)  and  (2.3)  with 

n(r  -1)  replaced  by  n(i-0  )  and 
n  n 


(3.1) 


THEOREM 


a  =  8  log[-(l-p  )/\p’(r  )]  B  =  -1/log  r  . 
n  n  nnn  n  n 


3.1.  Suppose  G^d/y^)  converges  weakly  to  Gq  with  Tq  =  1,  and  lim  ^  00 . 


Then  the  possible  limit  distributions  for  M  are  F  ,  0  <  c  <  °°.  The  M  has  limit 

n  c  —  —  n 

distribution  Fc  if  and  only  if  nd-p^)-^.  Appropriate  norming  constants  for 
the  convergence  (2.2)  are  given  by  (2.3)  and  (3.1). 

Proof.  This  follows  by  a  proof  paralleling  that  of  Theorem  2.1.  Here  the 
maximum  queue  length  in  a  busy  cycle  has  the  distribution 


•S'*.*  •  '  -  • 


■  ,•  ■  •  ■/  *  *  ■ 


■*.  .  ^ ■*,  ^ ^ ^ ™  ***  /*%<"  'i* »  * 


where  C  is  a  circle  in  the  z-plane  with  center  at  the  origin  an.  r  ius  less 

than  r  ;  see  Cohen  (1969).  Choose  circles  and  C  in  the  z-plane  with 
n  n  n 

centers  at  the  origin  and  radii  less  than  r^  and  greater  than  1,  respectively, 

such  that  i p  (z)  has  no  zeros  other  than  r  and  1  on  the  annulus  between  C^ 
n  n  n 

and  C  .  Then  by  the  residue  theorem 


(2iri) 


-1 


z  n\j>  (z)  ^dz  -  (27Ti) 

n  J 


z  "fy  (z)  ■'"dz  =  ^'(1)  ^  +  r  m^'(r  )  \ 
n  n  n  n  n 


,n 


n 

Using  this,  (3.2)  and  4>*(1)  ^=p  *"-l,  we  have 

n  n 


(3.3) 

where 


n[ 1  -  F  (a  +b  x) ]  =  {(£  -p  )/[n(l-p  )]  +  e  } 
n  n  n  n  n  n  n 


-1 


-m 

C  *  -r  n(l-p  )/4»'(r  ) 
n  n  n  n  n 


e 

n 


(27Tin) 


[  z  n4>  (z)  ^dz 


and  m  is  the  integer  part  of  a  +b  x.  Arguing  as  in  the  proof  of  Theorem  2.1, 
n  n  n 

one  can  show  that  n(l-Pn)  -*■  c  implies  (2.10)  and  hence  (2.2).  The  rest  of  the 
proof  is  the  same  as  before. 


REMARKS  3.2.  Analogues  of  Remarks  2.2  and  Corollary  2.4  apply  to  Theorem  3.1. 
Note  that  since  (3.3)  is  simpler  than  (2.13),  Theorem  3.1  does  not  need  an 
assumption  like  (2.5). 
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